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Abstract 

We examine the problem of the evaluation of both the propagator 
and of the partition function of a spinning particle in an external 
field at the classical as well as the quantum level, in connection with 
the asserted exactness of the saddle point approximation (SPA) for 
this problem. At the classical level we argue that exactness of the 
SPA stems from the fact that the dynamics (on the two-sphere S 2 ) 
of a classical spinning particle in a magnetic field is the reduction 
from R 4 to S 2 of a linear dynamical system on R 4 . At the quantum 
level, however, and within the path integral approach, the restriction, 
inherent to the use of the SPA, to regular paths clashes with the 
fact that no regulators are present in the action that enters the path 
integral. This is shown to lead to a prefactor for the path integral that 
is strictly divergent except in the classical limit. A critical comparison 
is made with the various approaches to the same problem that have 
been presented in the literature. The validity of a formula given in 
literature for the spin propagator is extended to the case of motion in 
an arbitrary magnetic field. 
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I Introduction. 



Since the early days of path integration, how to do a path integral for spinning 
particles was recognized as one of the major difficulties of the formalism. 
Schulman @(but see also 0) made a first attempt towards a formulation of 
a path integral for spinning particles, one which was however rather a related 
path integral, namely that for a spinning top. 

Much progress has been made since with the systematic use, initiated by 
Klauder ||, of the resolution of the identity associated with spin-coherent 
states |7j in the discretized-time (or time-sliced) approach to the path 
integral. Path integral quantization using coherent states will be discussed 
extensively in the sequel. Other versions of the path integral not making 
an explicit use of coherent states have been discussed instead in || [| (see 
also 13). 



Already at the classical level, a spin (classically a vector S of fixed magni- 
tude s) presents some peculiarities as a dynamical system. While the Hamil- 
tonian description is essentially straightforward if one assumes the Poisson 
brackets: 

{Si, Sj} = eijkSk (1) 
among the components of the spin vector, it has been pointed out by Bal- 



achandran et al. |]TT| that no global Lagrangian description can be given as 
long as one sticks to the natural configuration space of the spin, which is 
the compact two- sphere S 2 . The same authors have shown that a global La- 
grangian can be associated with a classical spin by lifting the description of 
the system from S 2 to the group manifold of SU(2), which is the three-sphere 
S 3 , as follows. Considering the usual spin-1/2 representation of SU(2), one 
can define a vector S in S 2 via the Hopf map, i.e.: 

S = S ■ a = sga z g~ l , (2) 



with a = (o"i, 02, 03) the Pauli matrices. Then it can be shown [LJJ that the 
(global) Lagrangian on TSU(2): 

L = isTi{a,g- l g) - ^Tr(SB) , (3) 

where fj, is the Bohr magneton and B = B ■ 0, yields the correct equations 
of motion for a classical spin in an external magnetic field B. The same 
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equations can be derived of course at the Hamiltonian level using the Poisson 
brackets ([!]) and the Hamiltonian: 



H = fiS-B. (4) 

As it is clear from the fact that the spin is recovered via the Hopf projec- 
tion (0), this approach introduces an extra, nondynamical U(l) gauge degree 
of freedom. Indeed, under: g — > g exp fry 03 /2] , (0) is invariant while L 



changes by a total time derivative (i.e. it is "weakly invariant" [IT]): 

L ->• L - 57 . (5) 
This is also evident if we parametrize SU(2) with the Euler angles as 

g — e -i^3/2 e -iea2/2 e -i-ya 3 /2 /g\ 

(0<#<7T,0<</>< 2tt, < 7 < 47t), which yields the (local) Lagrangian: 

L = s(0 cos # + 7) + fiS ■ B . (7) 



Although it can be shown [|12| that this extra gauge degree of freedom 
has interesting consequences for the path integral quantization of (|3|), and 
namely that it leads in a straightforward way to spin quantization (i.e.: 
2s/fr=integer), we would like to stress here the fact that (||) (or (0) for that 
matter) being linear in the time derivatives, a spinning particle is described 
as a constrained dynamical system. Canonical quantization requires then the 
use of Dirac's theory |T3] of constraints, and it has been shown by Balachan- 
dran et al. |ll|] that this does indeed yield the correct quantum mechanical 
description of the spin including spin quantization. As pointed out, e.g., in 
Ref. 0, it turns out that the (semiclassical) Bohr-Sommerfeld quantization 
is exact for a spin Hamiltonian. This is of course a strong indication that 
the saddle point, or stationary phase, approximation (henceforth referred to 
as SPA) to the path integral for, say, the propagator should be exact as well. 

Concerning the path integral approach and the legitimacy of evaluat- 
ing the path integral within the SPA, which is the main problem which the 
present paper addresses to, a Lagrangian of the form (0) poses another prob- 
lem, and namely that, the classical equations of motion being first-order in 
time, the saddle point problem becomes overdetermined, as the solutions of 
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the saddle point equations have to obey two boundary conditions and not a 
single one. Different ways to overcome this problem have been proposed in 
the literature, and notably by Klauder [|, Keski-Vakkuri et al. [[HJ] and 
Suzuki |fl5f| . 

This paper is devoted mainly to the discussion of two questions, namely 
to that of the legitimacy of some approximations that are currently made 
in the coherent state formulation of the path integral for spins and of the 
apparently surprising fact, that has also been widely discussed in the lit- 
erature, that the SPA to the path integral yields the exact result for both 
the propagator and/or the partition function of a quantum spin. While the 
SPA is almost trivially exact for quadratic Hamiltonians, here the apparently 
surprising fact is that it turns out to be exact also for the Hamiltonian 
and/or the Lagrangian (|7|) that are far from being quadratic. Exactness of 
the SPA implies of course that we discuss the applicability to the present 
problem of the Duistermaat-Heckman (DH) theorem which provides a rigor- 
ous framework for the discussion of the exactness of the SPA. 

The paper is organized as follows. In the Appendix we briefly review the 
DH theorem and the SPA. In Section II we apply the latter to the calculation 
of the partition function of a classical spin and show that the deep reason of 
the validity of such an approximation is that a classical spin, when viewed 
as a dynamical system, can be shown to result from the reduction of a linear 
dynamical system from R 4 to S 2 . In Section III we begin by discussing briefly 
the use of spin-coherent state path integral approach to the calculation of 
the propagator and/or of the partition function of a quantum spin. We show 
there that the approximation currently used in literature restricting the paths 
in the functional integral to be continuous leads to incorrect and diverging 
prefactors. Then we review briefly the approaches of Klauder and Suzuki 
to the same problem. We clarify the origin of some apparently mysterious 
terms that are added to the action in Klauder's approach. Then we analyze 
why the SPA yields the exact result in Suzuki's approach, by arguing that his 
main result stems simply from the exactness of the Ehrenfest theorem in the 
present context. In Section IV we discuss in more detail the coherent state 
approach by using the holomorphic representation for the latter, applying 
the complex SPA to the discrete version of the spin path integral. The final 
Section V is devoted to a general discussion and to some conclusions. 
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II The classical spin. 



Let us start our considerations from a classical (nonrelativistic) spin in a 
time-independent magnetic field B = (Bi, B 2 , B 3 ). 

Since the only degree of freedom for a spin is its direction, we describe 
it by means of a three dimensional vector (S%, S 2 , S3) E R 3 of fixed norm: 
Sf + S 2 + S3 = s 2 , where the three classical variables Sj (j = 1, 2, 3) satisfy 
the Poisson brackets ([!]). Thus the phase space for a classical spin is the 
two-dimensional manifold S 2 (s), equipped with the symplectic two-form 



Q = —e ijk SidSj A dS k . (8) 



The Hamiltonian that describes a classical spin in a magnetic field is given 
by (U). For simplicity we will set // = 1 in the following so that 

H = B-S. (9) 

The Hamiltonian vector field A associated to (H), and determined by the 
symplectic form Q, is given by 

A = ^kBjSk-^- , (10) 

so that the classical equations of motion, = dH, read 

Si = eijkBjSk ■ (11) 

Without any loss of generality, we will set B = (0, 0, B) (B > 0) from now 
on. In this case, the Hamiltonian (H) and the equations of motion (pT|) 
assume a very simple form if we use spherical coordinates S% = s sin 9 cos </>, 
S2 = ssin^sin^, S3 = scos#. They become respectively 

H = sB cos 9 (12) 

and 

cos99 = (13) 
sin 0(0- B) = . 
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The latter can be easily integrated and one sees that the classical orbits are 
circles parallel to the equator, the spin precessing about the magnetic field 
with a period r — 

In spherical coordinates, it is also very easy to compute the exact partition 
function for the Hamiltonian (191): 



Z d {p) = / Qe^ B . (14) 

Indeed one gets: 

ZJp) = s r d<p smed0e-P sBcose (15) 
Jo Jo 



2vr 

Jb 



e (3sB _ e -/3sB 



From (|15|) , one recognizes that the classical partition function can be written 
as the weighted sum of two terms, each given by the evaluation of e~^ H 
at the two critical points 6 = tt and 6 = of the Hamiltonian for which 
As=(oo±s) = 0- I n addition, it is not hard to check that the weights (±fg) 
are exactly the ones coming from the calculation of the contributions to the 
integral flT5] ) of the gaussian fluctuations around the stationary points of H. 

Everybody is familiar with such a result whenever dealing with (multidi- 
mensional) harmonic oscillators or in general with a quadratic Hamiltonian 
on the linear manifold R 2n . Even if the Hamiltonian @ in not of this kind, 
the SPA is exact as well. A spin in a magnetic field is in fact the simplest 



(nontrivial) application of the Duistermaat-Heckman theorem [|T^, 17 , which 
establishes under which conditions a phase-space integral, such as (0), can 
be evaluated exactly in the SPA. We refer to the Appendix for a review of 
the Duistermaat-Heckman theorem and for the proof of its applicability to 
the system of a spin in a magnetic field. Here we recall only that this result 
holds essentially for the two following geometrical reasons: 

1) the Hamiltonian ( |T2D is invariant under an [/(l)-action, given by rotations 
about the third-axis (i.e. the axis of the constant magnetic field); 

2) the associated Hamiltonian vector field, given by 

a = b HI +s '£) (i6 » 
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is proportional to the generator of this [/(l)-action. This is clear in spheri- 
cal coordinates, where A = B-^r. However while (^) defines it globally on 
S 2 (s), the sperical coordinate representation becomes singular at 9 = 0, tt. 



The Duistermaat-Heckman theorem gives some abstract mathematical 
conditions for the SPA to be exact. In the following we will show that for 
a spin in a magnetic field there is however a deeper reason why this holds: 
the dynamical system that describes a spin in magnetic field is the reduction 
of a bigger system which is described by a quadratic Hamiltonian on the 
linear manifold R 4 [18] . Such a situation has already been considered in the 
literature (see and references therein), even if not much in the context 
of the SPA, but mainly in the context of integrable systems, for which there 
exists the conjecture |19j that every integrable system is the reduction of 
a bigger linear dynamical system. To show what happens in the case of 
a spin, let us consider the linear manifold R 4 = {(xi, x 2 , x 3 , x 4 )} = C 2 = 
{(zi,z 2 ) : Zi = xi + ix 2 ,z 2 = x 3 + 2x4}, equipped with the symplectic two- 
form Cl = w-(dxi A dx 2 + dx 3 A cfe 4 ) = ^-^(dzl ^ ^Z\ + ^ z 2 ^ dz 2 ). 

There is a natural action of SU(2) on R 4 = C 2 , which is simply given by 
left multiplication: 



'a -(3*' 








azi — j3*z 2 


(3 a* 




. Z2 . 




(3 zi + a*z 2 



where 



a —(3 
13 a* 

mental representation. This action is symplectic and Hamiltonian 
its Lie algebra is spanned by the vector fields: 



, with \a\ 2 + \(3\ 2 = 1, is an element of SU(2) in the funda- 

, and 



Xi = xidi - x 3 d 2 + x 2 d 3 - xi(9 4 
X2 = —x 3 di - x A d 2 + x x d 3 + x 2 d A 
X3 = x 2 d\ - x\d 2 - Xid 3 + x 3 <9 4 . 



(17) 



Hence it is easy to prove that the linear vector field A = X^=i BjXj is a 
Hamiltonian vector field with a quadratic Hamiltonian, namely: 



H 



J2 B ifi ' 
i=i 



(18) 
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where 



/i = ^(^1^3 + x 2 x 4 ) = ^-Rez{z 2 
2s 2s 

h = y s ^ XiXa ~ X2X ^ = 2s lmz ^- Z2 

h = ^tf + 4-4-£) = ^{\*i\'-\*\ 3 ) • (is) 

Here we can identify R 4 with T*R 2 with canonical coordinates x%, X3 and 
momenta Pi = ff and P2 = ff- 

The SU(2)-action we are considering leaves the three-dimensional spheres 
S 3 (R) = {(xi,x 2 ,x 3 ,x 4 ) : x\ + x 2 + X3 + x\ = R 2 } = {{zi,z 2 ) : | ^1 1 2 -t- | ^2 1 2 = 
R 2 } invariant so that we can restrict the dynamics from the full R 4 to the 
submanifolds S 3 (R) on which the classical orbits lie. In the following we 
choose R = 2s and work on S 3 (2s). 

The functions ( |T9D and hence the Hamiltonian ([18]) have an additional 
symmetry, being invariant under the action of U(l) given by: 

(zi,z 2 )^(e i0 zi,e i9 z 2 ) 6e[0,2ir[. 

This allows us to project the Hamiltonian ( fL8|) from S 3 (2s) down to the 
two-dimensional manifold S 3 (2s)/XJ(l), which is homeomorphic to the two- 
dimensional sphere. Indeed, the three functions (x\, x 2 , x 3 , x 4 ) 1— > Sj = 
fj(xi,x 2 , x 3 , X4) (j = 1,2,3) given in (|19D are the components of the pro- 
jection map from the three-sphere S 3 (2s) to the two-sphere S 2 (s) of the 



Hopf bundle U(l)-> S 3 -> S 2 |l 



On S (s) the Hamiltonian (0) becomes simply H = Ylj=i BjSj an d there- 
fore coincides with the Hamiltonian (H) for a spin in a magnetic field. This 
shows that the latter is the reduction of a quadratic Hamiltonian defined on 
a linear manifold [J 



1 To prove this rigorously, we should show also that the restriction to S 3 (2s) of the 
symplectic two- form Cl is the pull-back of the symplectic two-form fl we defined on S 2 (s) . 
The proof of this statement is straightforward, so we will omit details here. 
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Ill Coherent state path integrals for spin. 



Let us consider now the quantum mechanics of a spinning particle, described 
by the Hamiltonian 

H = B-S, (20) 

where the spin operators Sj (j = 1,2,3) satisfy the usual commuation rela- 
tions {h — 1): 

[Si,Sj] = ieijkSk ■ (21) 

We have already mentioned in the Introduction that the semiclassical Bohr- 
Sommerfeld quantization turns out to be exact for this problem and this 
seems to suggest that the SPA to the path integral for the propagator and/or 
the partition function could be exact as well. In addition, for such a simple 
problem, one can evaluate the partition function exactly. Its expression 

p f3Bs p -/3Bs 

Z((3) = Tr{e-^} = T -^s + ( 22 ) 

can be thought of as the sum of two terms, each corresponding to one of the 
two poles of the sphere (S 3 = ±s), similarly to what happens in the classical 
case. 

In this Section we will briefly review the existing literature on the evalu- 
ation of the partition function for a quantum spin in a magnetic field and in 
particular on the validity of the SPA applied to this problem. 

Let's begin by considering a set \l) of generalized coherent states 0, [7| 
labeled by one or more continuous variables that we denote collectively as I. 
The |Z)'s will be assumed to be normalized. They are however overcomplete, 
because although there is a resolution of the identity associated with them: 

1 = J dl\l)(l\ (23) 

with ll dl n a suitable measure, in general they fail to be an orthonormal set: 
(/|/')^0. 

Here we are interested in the group SU(2), which is generated by the spin- 
operator algebra (|2T|). Coherent states [0 for a spin s (2s being an integer) 
can be constructed as 

\e,<j>) = e -^ 3 e-^ 2 |0) , (24) 
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where |0) denotes the highest-weight state of the spin-s representation of 
SU(2) (53 1 0) = sjO)) and 9, are the two angular coordinates parametrizing 
S 2 : < 9 < vr; < <f) < 2tt. It is well known §, [7| that: 



10, 



+ sin- sin -e'^'-^ 2 
2 2 



2> 



(25) 



and that the resolution of the identity associated with spin-coherent states 
is: 



2s + 1 r 2 * r 



An Jo Jo 



' if \ ) 



sm 



(26) 



Going back to the general case, let H be the Hamiltonian of a quantum 
system. The matrix element of the propagator 



K(l F ,l I ,T) = (l F \e- iHT \l I ) 
can be represented as the following path integral: 



N N 

K(l F ,l!,T) = lim / n(4) II<Wi|e- ie ^ fc ) , 

"' fe=l fc=0 



(27) 



(28) 



where e = T/N and |/ ) = |^/), \In+i) = 1(f)- Eq. (p8|) can be rewritten as: 



iA(l k+1 ,l k ) 



where 
and 



r N N 

K(l F ,l r ,T) = lim / U(dl k )l[e 

J k=l k=0 

A(l k+ x,l k ) = —i^-ih+ilh) — e H(l k+ i,l k ) 

TT ( J , \ (h+i\H\l k ) 
ti KJk+l-, L k) 



(29) 
(30) 

{h+l\l>k) ^ ^ 

It is usually assumed that, for e small, is so close to \l k ) that one is 

allowed to expand the former around the latter to leading order in e. This 
leads to 

H(l k+1 ,l k ) ~ H(l k ) = (l k \H\l k ) (32) 

and to 

ln(l k+1 \l k ) ~ (5l k \l k ) , (33) 
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where \8lk) = — Hence the continuum version of the path integral 

reads 

rT 



Vlexp 



i I dt 

o 



%{l\l)-H{l) 



(34) 



where [/) = d\l)/dt. This procedure is justified when the Hamiltonian con- 
tains an explicit kinetic term that can act cts db rc gulator concentrating the 
functional measure on continuous paths. This is the case of Wiener inte- 



gral, where [22], |23| the measure is concentrated on paths q(t) satisfying the 
Lipschitz condition: \q{t + e) — q(t)\ = 0(e l l 2 ). The same holds true for the 
Feynman path integrals for massive particles in not too singular potentials . 

In the case of spin, the Hamiltonian fl2~0|) contains no regulators and the 
straightforward application of ( |3"2"D and ([£J) is questionable. We will see 
indeed that it may lead to serious problems. Let's be more specific. Taking 
again the magnetic field along the positive z-axis we have, in terms of spin- 
coherent states (|2l 



(l\H\l) -> {6,(j>\H\e,(j>) = sBcos9 



and 



i(l\i)^i(9,0\±\9,, 



s cos( 



So we obtain the path integral: 
(0 F ,4) F \e- illB ^ T \e I ,4> I ) = [vVlexp 



s cos 9(f) — sfiB cos 6 



dt 



(35) 
(36) 

, (37) 



where T>Vt = sm9T>9T)(f). From (p7[) we can obtain the canonical partition 
function 2(/3) by setting T = — z/3, \9f,4>f) — \@i,(j>i) an d by tracing over 
the angles. Following the standard procedure and Ref. 0, we have evalu- 
ated Z{(3) by time-slicing the path integral and using the resolution of the 
identity (p6|) at each intermediate (Euclidean) time. Denoting with Z N (j3) 
the result for time slices, we find for it without any further approximation 



the expression [12 



Z N {(3) 



1 + 



1 



N 



(38) 
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where Z(f3) in the right hand side is the exact partition function (p2|)p|. 

As one can easily notice by simple inspection, in the limit iV — > oo, Zn 
coincides with the exact partition function up to a diverging prefactor. Such 
divergency disappears only in the classical limit, namely h — > 0, s — > oo with 
hs = const. On the other hand the only approximations we used are those 
given by eqs. fl3~2"|) and mainly (j33f ) which again are exact in the classical 
limit. In the following of the paper we shall come back to this point. 



We discuss now briefly two different approaches g, [T5| that lead to the 
conclusion that in the case of a spinning particle in a magnetic field the SPA 
to the path integral for the propagator does indeed yield the exact result. 
This holds again up to a normalization factor not taken into account in[||, |T5| ], 
since the calculation is done in the continuum. 

Klauder |J[ has proposed a modified form for the action that appears in 
the path integral fl3"3|), redefining the latter as: 





fT 


/ VI exp 


i dt 




Jo 



i(l\l) + -ie(i\(l-\l)(l\)\l)-H(l) 



39) 

The prescription here is that the limit for e — > should be taken after having 
evaluated the path integral in the SPA. For the case of a spin in a magnetic 



field, (|3~9p becomes: 



e^O 







j VQexp 


['/*[ 

Jo . 



s cos 9(f) + -ise(0 2 + sin 2 9(j) 2 ) — sfiB cos 9 



(40) 
dt 



and the action in ( ]4C|) becomes essentially that of a particle of charge s and 
mass m = l/2se moving on the two-sphere S 2 , coupled both with a magnetic 
monopole of unit strength located at the centre of the sphere and with a 
constant electric- type field directed along the z-axis. In this context, Dirac's 



quantization condition |24j becomes identical with the spin quantization con- 
dition, i.e.: 2s =integer. 



2 In order to include paths turning around north pole any number of times we need 
to extend, at every time slice, the domain of integration of the variable <fi from [0, 2ir] 



to [— oo, +oo]. In fact without following such procedure one would get, as in Ref. 
completely wrong partition function. 



a 
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The classical equations of motion that can be derived from the action 
in fllOD are now second-order in time for e 7^ 0, and hence are not plagued by 
the already mentioned problem of overdetermination. Klauder has proved 
that they can be solved explicitly and that the resulting SPA to the path 
integral (or, better, what he calls the "dominant SPA" , namely approximat- 
ing the path integral with e iScl , without any prefactor originating from the 
integration of gaussian fluctuations) is indeed exact. Actually, only the free- 
spin case B = has been considered in Ref. ||, but the extension to B 7^ 
is straightforward [0. It should be noted, incidentally, that the term that 
Klauder has added to the action is in the form of a kinetic energy term, thus 
providing the required regulator justifying the assumption of continuously- 
varying paths. 

The origin of the additional kinetic-type term, that looks somewhat mys- 
terious in Klauder's original paper M , can be clarified as follows. If we push 
the expansion of the logarithm in (|33| ) one step further beyond first order we 
obtain, with \5l) ~ e\l) + e 2 /2|Z>: 



e 2 •■ e 2 • 

\n(l k+1 \l k ) ~ e(l k \l k ) + —{h\h) ~ -^{{k\k)) 2 



leading, in the continuum limit to 



lim exp 



/ / dt 

u 



i(i\i) + -ie(-(i\i) + (i\i)(i\i) 



(41) 



(42) 



where we have used — A final integration by parts of the term 

containing the second derivative yields precisely Klauder's additional term. 



Suzuki |T§] has adopted quite a different approach. Introducing two ad- 
ditional resolutions of the identity, the propagator 



where 



K(l F ,lj,T) 
A FI (l) = 



l(T)=h 
l(0)=h 



Vle iApi{l) 



dt 



i(l\l)-H(l) 



is rewritten in Ref. [13] in the form: 

K(l F ,l T ,T) = f I ' dl f dk{lF\lf){l f \e- lflT \k){h\li) 



(43) 
(44) 
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// 



A(T)=lt 

dl f dk(l F \l f )(k\li} / Vie"*® . (45) 

Jl(0)=k 



In general the parameter(s) I labeling a coherent state can be defined in terms 
of the expectation values of a suitable set of operators: position and momen- 
tum for the free particle and the harmonic oscillator, spin components in the 
case of spins. Taking then the latter as Cauchy data for the canonical equa- 
tions of motion one can determine how they evolve classically in time, thus 
determining a "classical" coherent state \l(t)) as the coherent state labeled 
by l(t), the evoluted at time t of any given initial parameter set I. Then, 
Suzuki has argued that the SPA to the last integral leads to 



r l(T)=lr 

/ Vle iA f M ~ 6(lf - k{T))e lA ^ , (46) 



where k(T) is the evoluted at time T of and that the final result of the 
application of the SPA is the semiclassical propagator 

K sc {l F ,l T ,T) = J dlS F UT)){k\h)e iA ^ , (47) 

where Au is the classical action evaluated along the classical path leading 



from 1(0) = U to l(T). Once again, the semiclassical propagator (j47|) yields 
the exact result for the propagator of a spinning particle. 

We would like to show here that the possible exactness of Eq. (f4~7D has 
very little to do with the path integral formalism itself, and that it follows 
rather from a single assumption, one that amounts basically to assuming 
that the Ehrenfest theorem be applicable in the present case. Let I be the 
set of operators whose expectation values label a given coherent state. To 
each coherent state \l) we can associate two different time dependent states: 

- the "classically" time evoluted state \l(t)) defined above and therefore such 
that 

(l(t)\l\l(t))=l cl (t), (48) 
l c iit) being the classical trajectory, 

- the quantum time evoluted state \l,t) obtained through the application of 
the full quantum evolution operator: 

\l,t) = e-^\l) . (49) 
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We require now the expectation values on the quantum evoluted state to 
evolve in time according to the classical equations of motion: 



(l,t\l\l,t) = l d (t) . (50) 

In other words we are assuming, as anticipated, the validity of the Ehrenfest 
theorem. A sufficient condition for this to happen is that the Heisenberg 
equations of motion be linear equations. This is certainly true for quadratic 
Hamiltonians with conventional canonical coordinates and momenta and the 
standard Poisson brackets among them. For spins, the Hamiltonian (|J) is 
not of the quadratic type, but the Poisson brackets (p}) lead nonetheless to 
linear equations of the motion and hence to the validity of the Ehrenfest 
theorem. If we now assume that the set of operators / act irreducibly on the 
Hilbert space of states (which is true, e.g., for harmonic oscillators and for 
spin systems) the two states \l,t) and \l(t)) will differ at most by a phase, 

\l,t) = e^\l(t)) . (51) 
Differentiating with respect to time we obtain 

m\j t \m) = ^x+(i,t\j t \i,t) (52) 

and hence 

(l\i) d = -ix-i{l,t\H\l,t) . (53) 

If, as it is presently the case, the Hamiltonian does not contain time deriva- 
tives, we obtain furthermore: 

X = i(l\i) d - (l(t)\H\l(t)) = i(l\i) d - H{l{t)) (54) 

and eventually, integrating the last equation and up to an irrelevant constant 
phase: 

x (t)= Hiiiiiu-Hiijlde . (55) 

JO 1 J 

The r.h.s. of ( [55"D is exactly the classical action, hence x{t) = A c i(t). Then 
Suzuki's result follows at once from: 
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K{l F ,k,T) = (l F \e- iHT \k) 

= J dk{l F \e- lflT \k){k\li) 

= J dl l (l F HT)){k\l I )e^ 

= JdhihMTMilkW^ , (56) 

which is the desired result. 

Let's analyze fl56|) in the specific case of a spinning particle. As already 
pointed out the Hamiltonian (|J) leads to linear Heisenberg equations of mo- 
tion and hence to the validity of formula (0). It seems to us that this is the 
ultimate reason for the apparently surprising fact that SPA to path integrals 
and/or semiclassical quantization leads to the exact result for the quantum 
propagator of a spinning particle. We would like to stress here that this is 
also true for the motion in a magnetic field with an arbitrary dependence on 
time. Formula (|56"D reduces therefore the quantum problem of the calculation 
of the spin propagator in an arbitrary magnetic field to the solution of the 
classical equation of motion. 

The expression for the propagator can of course be worked out explicitly 
only when the classical equation of motion can be solved analytically. With 
a constant magnetic field along the z-axis, formula (|47|) reads: 

K(8 F ,4> F ,9i,4>r,T) = J ^(^,0 F |ft t (T),0 4 (T))(^,0^ / ,0 / )e^ . 

(57) 

Since the solutions of the classical equations of motion are 0i(T) = 9i\ 4>%{T) = 
(pi + fiBT, while A c i = 0, the above expression becomes: 

K{p F ,<t> F ,6i,<t>f,T) = J da<0j^|0i,&+/iST)(0 i > i |0 J > J ) . (58) 

The integral in ( |58D can be done easily by noting that 
(9 F , 4> F \9i, 4>i + fiBT) = (9 F ,(f) F — fiBT\9i, (pi) and using the resolution of 
the identity. We obtain eventually 

K(9 F , (P F} 9j, (Pr, T) = (0 F , 4>f - \iBT\Q u fa) , (59) 

which indeed gives the exact propagator, upon using formula (|25|) for the 
overlap of two coherent states. 
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IV The complex saddle point. 



We devote this Section to another approach to path integral for spin. It has 



been used for example in |L4[ and in ||25|| , where however only continuous 
and hence formal expressions for the action and the path integral have been 
studied. 

Here we will study the path integral discrete version for the propagator 
fl34|) . The method is based on a different choice of the spin-coherent states. 
Instead of (P^D, we consider the following states, labeled by a complex pa- 
rameter fi: 

e fiS~ 

l"> = (TTwF |0> ' (60) 

where we have set S' ± = Si ± i£>2- To be more precise, \i parametrizes 
the homogeneous space SU(2)/U(1) = S 2 by means of the stereographic 
projection of S 2 onto R 2 . With respect to the spherical coordinates 9 and </> 
used previously, one has \i = tan|e 1 ^ and it is also not difficult to check that 
the state flBOD coincides up to a phase with the state \9, 0) defined in 

To write down an explicit expression for the path integral in terms of the 
coherent states (|6"0|), we need the matrix elements 

(1 + A» 2 * 

(%) = (i + M»)-(i + |A|g)« (61) 

1 + A*/x 

W) = , , I | 2 (/V - A» 

and the explicit formula for the resolution of the identity: 

2s + 1 f d 2 u . . . . „ 

—JiTTWr^^ 1 (62) 

where d 2 fi = d(Re/x)d(Im/x). After some algebra, one finds that the path 
integral representation of the propagator is given by 

K(^r,T) = J MQH/i/ (T^pj^. (63) 
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where the action for a spin in a magnetic field is 




2 



(64) 



1 + |/i 

Again the SPA equations that are derived by extremizing fl6"4|) : 

{ r= -V < 65 > 

are first-order, so that they do not admit of a solution for boundary conditions 
of the form: 

f Mo) = m m 

with arbitrary fij and ftp. 

This problem does however admit of solutions |26| if we enlarge our 



variable space from C to C 2 and look for saddle point solutions for which 
fj,(t) and /U*(t) are not complex-conjugate one of the other. In other words, 
we have to consider /z and /i* as independent complex variables and look for 
solutions of (|65D satisfying the boundary conditions 

1 V (T) =H* F 

These so called complex saddle point solutions are easily found to be 

elBt ^ 

= e lB ^fi F m 
and, to complete the SPA, one has further to evaluate the contributions of 




the gaussian fluctuations around the solutions (^S 



This has been done by the authors of [III, E31, who have concluded that 



one can obtain the exact propagator in this way. They have however worked 
with the continuous expression ( p^ ) and therefore they have performed only 
formal calculations. Funahashi et al. |27 have considered the discrete version 



of the integral but they have applied SPA to the calculation of the partition 
function only. This problem, as pointed out by the same authors, is simpler 
than the calculation of the propagator since the overdetermination problem 
can be overcome without enlarging the variable space from C to C 2 . 
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Our aim in this Section is to examine whether the SPA to the path integral 
(63) is exact by performing the calculation in the discrete version of the path 
integral. The latter is given by fl2T3|), where now the generic index / stands 
for the complex parameter fi. By using ( |BTD and rearranging the terms, one 
eventually gets: 

-tf/n 1 ^-^ ,69) 

with 

*=r^r J., po) 



TT J (1 + |/1 F | 2 ) S (1 + \fIj\ 2 Y 

whereas the discretized action A is given by 



7V 



A = E [^iJ-i + H M-i] + 2s ln (! + VfUn) (71) 



A-,_ 1 = 2 a ln 1 + ^ J '" 1 (72) 

= ~^ se i 7 ,1, • ( 73 ) 

The expression for Ajj-i given in ([72] ) comes from the exponentiation 
of the overlaps (j = l,---,iV). Let us remark that we have no£ 

expanded the overlaps to first order in the difference — l/^-i), as usually 
done to recover the continuum version (see (|33D ). As already pointed such an 
expansion would be correct only in the presence of a regulating term in the 
action. Thus we will work only with (|72|) and we will proceed now to evaluate 
the SPA to the multidimensional integral (|69|), according to the formula 

N-l / i \ JV-1 

where Mn is the 2(N — l)x2(N — 1) matrix of quadratic fluctuations around 
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the classical solutions of the discretized action ([H|): 



A = A d + (S/il, 6/j.i, 5n* N _ x , 5li n _i)M n 



( $Vi \ 
6(4 

V W-i J 



(75) 



The subscript "c/" means that the function (the action in this case) has to be 
evaluated on the classical solution, which satisfies the saddle point equations: 



Mj-Mj-i 



-ieB- 



"j+l 



+ieB 1 /r 1 

l+KA»j-i 



(76) 



In ([76|), as well as in ([H]), we have explicitely written \fij\ 2 as Lijii* to stress 
the fact that in the search of the saddle point solutions we have to treat [Aj and 
/i* as independent complex variables. Indeed, exactly as in the continuum, 
the equations flT6|) are incompatible with the boundary conditions 



yUo = Hi 



(77) 



unless Lij and //* are treated as independent variables. 

We know that the classical solution has to fulfill the property Lij+i\ c i = 
f^jlci + 0(e), so that we can approximate the denominators in the right hand 
side of ([76|), to get 



1+fj.jH 



+ Ol 



0(e 2 ) 



(78) 



To order e, we find for the solutions satisfying the boundary conditions (ff"7|) 
and for the classical action the expressions: 



[1 + ieBYm 
ll + ieB) N ^/i* F 



A 



cl 



-isBT 



(79) 
(80) 
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To complete the calculation of the right hand side of ( |74|) we have to 
compute the determinant of the gaussian fluctuation matrix Mjy. As for the 
solutions of the classical equations flTff), we evaluate the matrix elements to 
the order 0(e). Neglecting therefore all the terms of at least order 0(e 2 ) we 
find that the only non-zero matrix elements are 



d 2 A 



-2s 



d 2 A 



el 



1 + fifij 

d 2 A 



-2s 



cl 



[1 + (1 + ieB) N fi* F fi T ) 2 ' 
d 2 A 



2s 

1 + 



1 



1 + tflJlj 



+ ieB 



1 



(1 + ^ +1 Hjf 



2s(l + ieB) 



{\ + {l + ieB) N ^ F ) 2 



Defining 



A, 









ft; 



B; 




h 



hi 




the matrix is given by 

/ Ax B x 



M, 



N 



C\ A2 B2 



c 



N-3 



An-2 B 



N-2 



Cn-2 A n _i J 



Thus we have 

detMiv 



b N-l 



-2s 



2(JV-1) 



l + (l + ieB) N fijn* F ) 2 _ 
which can be finally inserted, together with (0), in (|74l) , yielding: 



K(fJ,* F ,fMi,T) t 



N-l 



-isBT 1 



1 + mn* F {l + iBT/N) 



N\2s 



;i + M 2 Hi + M 



(81) 



(82) 



(83) 



(84) 



(85) 



(86) 
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This is again the expected result up to the divergent normalization factor 
(1 + l/2s) N ~\ Indeed since lim^o^r + iBT/N) N = e iBT ', but for the 
prefactor, we obtain 

Kfr F , ^ T) sc = (1 + M2)s(1 + |/i/|2)s , (87) 
which coincides with the exact propagator. 



V Conclusions. 

In this paper we have examined the problem of a spin in a magnetic field. We 
have seen that the SPA applied to the calculation of the classical partition 
function yields the correct result. In quantum mechanics the situation is 
more complicated. We have reexamined the different approaches that have 
been used in literature to prove the exactness of SPA in the calculation of the 
quantum propagator and partition function. All these methods, in particular 
those proposed by Klauder || and by Keski-Vakkuri et at |14|1 , make use of 



the continuum expression of the path integral and hence reproduce the correct 
result only formally. 

To test the validity of the SPA for a spin system we decided to work with 
the very definition of the path integral, namely with its discrete version. We 
have written the discrete path integral for the partition function and the 
propagator in terms of spin coherent states and then considered two differ- 
ent kinds of approximations. In Section III we have expanded the overlaps 
(lk\h-i) t° first order in e = T/N and then performed an exact integra- 
tion. On the contrary, in Section IV we have made no expansion for the 
overlaps and applied instead the complex saddle point method to evaluate 
the SPA of the path integral. In both cases, we have found that the exact 
result is reproduced correctly only up to an infinite normalization factor, 
limjv^l + l/2s) N , which goes to 1 in the classical limit s — > oo, provided 
the latter is performed first. Notice that this fact is not a drawback just of 
the spin Hamiltonian. Indeed one can repeat easily the calculation for H = 
and get the same infinite constant. 

All this seems to suggest that, for those Hamiltonians that do not contain 
a regulating term (such as a kinetic part), any approximation that restricts 
the class of quantum paths in phase space, by imposing some regularity 
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conditions on them, yields a wrong result. In our case this shows up in the 
appearance of an infinite prefactor. 

We would like to conclude by commenting on the paper by Funahashi et 



al. who have been able to reproduce the exact partition function (with 
the correct prefactor) by performing an SPA in the discrete path integral, 
written again in terms of spin coherent states. They have obtained this result 
by taking into account also the gaussian fluctuations coming from the factor 
1/(1 + |/x| 2 ) 2 = sin 9/4 appearing in the integration measure. We do not want 
to explain this technique which is described in detail in [27|. Here we notice 



only that the inclusion of fluctuations coming from the measure induces a 
shift in the multiplicative factor appearing in front of the action from 2s to 
2s + 1. Effectively, we can say that such a method amounts to choosing 2s + 1 
as parameter for a semiclassical expansion. If so, we should not evaluate the 
path integral following (|73|), but according to: 



N-l ( 

A/-n 

3=1 \ 



A 



71 



N- 



-l e (2s+l)A ci 



(1+^)2 ( detM 



cl 



IN 



1/2 



where A = A/2s and Mjv is the matrix of gaussian fluctuations of (2s + 1)A 
around the classical solution, so that 



det M, 



N 



<W-1 



-2(8 + 1) 



[l + (l + ieB) N fi lf i F ) 



2(JV-1) 



(89) 



It is exactly the factor 2s + 1 in ( |39"D which cancels the same factor in A/", 
yielding the correct propagator. 

This is quite a remarkable result. In our opinion, however, the derivation 



presented in [g7j needs some clarification. To us it seems to be inconsistent to 
include gaussian fluctuations of the measure factor in the calculation of the 
SPA, without considering also its contributions to the saddle point equations. 
In [^] this method works only because the first derivatives of the measure 
factor do vanish at the classical solution, which in this case correspond to 
(j, = or /j, = oo (6 = 0, 7i in angular coordinates). But this would not be the 
case in slightly more complicated situations, for example for the calculation 
of the propagator or when considering complex saddle point solutions. 
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A The Duistermaat-Heckman theorem. 



Let us consider an oscillatory integral of the kind 

'«>-(£)" Jm"* (A1) 

over a (2n)-dimensional manifold M with volume form a. If M is a Rie- 
mannian manifold, under rather mild hypotheses, namely that the function 
/ be a Morse function, i.e. that the Hessian matrix of / be non-singular 
at all critical points of / (det Hessp(/) 7^ if V/(P) = 0), it is possible to 
show [E3] that for large values of the parameter t, one has 



I(t)=Y J c P e u ^ + 0(t- l ) 1 (A2) 
p 

where the sum ranges over all critical points of / and the coefficients are 
given in terms of the determinant of the gaussian fluctuations of / around 
the critical points: 



cp = exp 



7T 

i-sgnHess P (/) 



[detHess P (/)p . (A3) 



Here the signature sgnA of a symmetric real-valued nonsingular matrix A 
is defined as the number of its positive eigenvalues minus the number of its 
negative eigenvalues. 

Everybody is familiar with the elementary result that the remainder term 
0(t _1 ) vanishes identically if M is the linear manifold R n with volume form 
a = dx\ ■ ■ ■ dx2n and the function / is a quadratic form: / = \Qx ■ x — £ 



x. 



Q being any symmetric real-valued (2n)-dimensional non-singular matrix. In 
this case the only critical point of / is xq = Q _1 £ and Hess£ (/) = Q, so that 
l2]) with 0(t _1 ) = gives simply the formula for a gaussian integral. 
The theorem of Duistermaat-Heckman |TB[ and its generalization due to 



Berline and Vergne |17| establish under which conditions an integral of the 
kind (|Aip can be exactly evaluated in the SPA, i.e. when 0(t _1 ) = 0. 

Let M be a compact (2n)-dimensional manifold with symplectic two-form 
Q and suppose M is acted upon by a compact Lie group G, whose action 
is symplectic and Hamiltonian [p8|| . Let us denote with Xv the fundamental 
vector field on M generated by the action of the element 77 in the Lie algebra 
G of G and with f v the associated Hamiltonian function (i Xv Q = df v , where 
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% denotes the contraction). Then, if Xv * s non-degenerate, i.e. if it is zero 



only at the fixed points of G, the following results hold |ge 

1) f v is a Morse function; 

2) 

t \ n r ft' 
,27T/ Jm n! 



(A4) 



where the the sum ranges over the critical points of /, i.e. over the points P 
such that Xri(P) — and the coefficients cp are given by 



c P 



A1A2 • • • A n 



(A5) 



the Aj's being the coefficients appearing in the matrix Lp of the derivatives 

-g^ 1 J , which in a 
suitable positively oriented basis can always be written as | 



Ai 
Ai 



A 2 
-A 2 





-A r 







From a physical point of view, (|A4j) can be applied to classical statistical 
mechanics for the computation of the partition function: 



Z = [ ^-e-? H 
jm n\ 



(A6) 



In this case the vector field \v ^ s gi yen by the Hamiltonian vector field Ah, 
where now H plays the role of the function f v . The assumptions of the 
Duistermaat-Heckman theorem require Ah to be non-degenerate and to be 
the fundamental vector field associated with an element 77 of the Lie algebra 
G of a compact Lie group G acting symplectically on the manifold M. If 



If M is equipped with a Riemannian metric, formula flAq) is seen to coincide with 



(A3). 
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these two conditions holds, we can apply formula ( |A4| ) to conclude that (we 
have set (3 = —it): 

m = {j^J^c P e-^, (A7) 

where 

i) the sum ranges over the stationary points of the Hamiltonian, A#(P) = 0; 

ii) the coefficients cp are given by ( |A5|) in terms of the Ay's, the latter being 



the coefficients of the matrix [Xp]ij = ( ) , written in a suitable positive 
oriented basis. 



p 



Let us go back now to the problem of a spin in a magnetic field. We 
want to show that the Duistermaat-Heckman theorem can be applied to this 
problem, so that the exact partition function can be calculated by means of 
formula (|A7|). The Hamiltonian vector field associated to the Hamiltonian 
©is 

a = b H^ +Si £) (A8) 

and it is easy to recognize that it is proportional to the generator Xv = 
of the rotations about the third axis (i.e. the axis along the constant 
magnetic field). Thus the Lie group that acts symplectically on the phase 
space manifold S 2 (s) is simply given by U{1) in this case. 

To apply (|A7|) we have first of all to find the critical points of the Hamil- 
tonian, which are given by the North and the South poles of the sphere: 

P± = (0,0,±s) (A9) 

and then to compute the coeffcients cp ± according to ( |A5| ). In the tangent 
space of P + and P_ we choose to work with the positevely oriented basis 
(dfi ' 5^) anc ^ (afj ' a§r) res P ec tively. With respect to these bases 

L P± (A) = T B _° x I , (A10) 

so that 

cp ± = ■ (All) 
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We can finally compute the partition function in the SPA as 

Z = || (e^ - e-^) , (A12) 

which, in agreement with the Duistermaat-Heckman theorem, coincides with 
the exact partition function (|15|) for a spin in a magnetic field. 
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